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摘 要 : 针对 传统 回溯 搜索 优化 算法 存在 收 你 速度 


验 设计 的 改进 算法 。 首 先 将 正 交 实验 设计 方法 引入 算法 的 交 又 算 子 中 ,得 出 具有 代表 性 的 优质 子 代 个 体 ; 


胞 自动 机 邻居 模型 的 基础 上 ， 对 个 体 展开 领域 内 多 父 代 正 交 
与 交叉 的 种 群 引入 动态 优秀 个 体 比 例 权重 进行 选择 更 新 ， 并 采用 新 的 动态 变异 方程 ， 
索 能 力 。 通 过 对 12 个 标准 测试 函数 进行 仿真 实验 ， 并 与 其 他 6 种 表现 良好 的 算法 进 
法 在 收敛 速度 以 及 寻 优 精度 方面 都 具有 明显 优势 。 
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交叉 操作 ， 提 高 算法 的 开采 能 力 和 搜索 效率 ; 
平衡 算法 的 全 局 搜索 和 局 部 所 


行 比较 ， 结 果 表 明 ， 
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慢 、 搜 索 精 度 不 高 等 问题 ， 提 出 了 一 种 基于 元 胞 自动 机 和 正 交 实 


然后 在 元 
最 后 对 参 


改进 的 算 


Abstract: According to the Slow convergence Speed and low searching precision of traditional backtracking search 


optimization algorithm，this paper proposed an improved algorithm based on cellular automaton and orthogonal 


experimental design. Firstly, the algorithm introduced orthogonal experimental design method into the crossover operator to 


obtain representative high-quality offspring individuals; Then based on the neighbor model of cellular automaton, the 


method carried out the orthogonal crossover operation of multiple parents in the domain for individuals, which was 


beneficial to improve the mining capacity and search efficiency of the algorithm; Finally, in order to balance the global 


Searching and local searching ability of the algorithm, the method introduced the dynamic proportional weight of excellent 


individuals into the cross-population to select and update them，with adopted a new dynamic variation equation. The 


Simulation experiments Selected 12 standard test functions and compared with 6 other well-behaved algorithms, the results 


show that the improved algorithm has obvious advantages in convergence Speed and optimization accuracy. 


Key words: backtracking search optimization algorithm; cellular automaton; orthogonal experimental design; multiple 


parent orthogonal crossover; dynamic variation equation 


智能 优化 算法 是 受 人 类 智能 、 生 物 群 


象 规律 的 启发 而 提出 的 随机 搜索 算法 。 因 


并 


通用 性 和 易 操作 性 等 特点 ， 近 年 来 在 生产 


rs| 


复杂 的 约束 优化 问题 中 发 挥 了 重要 作用 
但 是 ， 由 于 BSA 算法 变异 算 子 中 


了 站 上 


功能 用 于 引导 搜索 ， 因 而 具有 较 强 的 全 局 搜索 能 力 ， 在 求解 


的 历史 种 群 引导 学 


于 


以 及 交叉 算 子 中 两 种 交叉 策略 随机 调用 都 


~ 


有 一 定 的 


目 


数据 挖掘 等 众多 工程 领域 得 0 泛 的 应 


能 优化 算法 主要 有 遗传 算法 ! 


] (genetic algorithm， GA)、 六 
子 群 算法 和 (particle swarm optimization，PSO)、 差 分 进化 算 
法 BCdifferential evolutionary, DE)、 人 工蜂 群 算法 外 (Cartifici 


bee colony，ABC) 等 ， 以 及 一 些 相 应 的 改进 算 光 
Liang 等 人 提出 的 粒子 群 改进 算法 CLPSO， 夏 学 文 等 人 [9 
改进 的 RLPSO; Wang 等 人 中 提出 的 OXDE， 周 晓 根 等 
提出 的 DELLU; 以 及 Gao 等 人 外 基于 蜂 群 算 光 
MABC， 周 新 宇 等 人 MM 改进 的 ABC-OED 等 。 
溯 搜 索 优 化 算法 0 (backtracking search algorithm ， 
BSA) 是 Civicioglu 在 2013 年 提出 的 一 种 新 型 
法 。 该 算法 由 于 结构 框架 简单 ， 控制 参数 较 少 和 独特 的 记忆 
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性 。 使 得 算法 拥有 较 好 的 探索 能 力 ， 但 开发 外 


E 力 较为 欠缺 ， 


在 求解 多 峰值 问题 时 算法 的 局 部 搜索 能 力 较 弱 ; 


optimization algorithm based on cellular 


CA-DBSA)， 分 别 从 交叉 算 子 和 变异 算 子 入 
中 引入 了 正 交 实验 设计 (orthogonal experimental design， 


改 敛 0 区 。 针 对 这 一 问题 ， 本 文 提 出 一 种 基 
动态 回溯 搜索 优化 算法 (dynamic backtracking search 


automata 


容易 早熟 
于 元 胞 自动 机 的 


9 


， 在 交叉 算 子 


OED) 方法 ， 并 在 元 胞 自动 机 (cellular automaton，CA) 的 


= 


邻居 模型 基础 上 进行 并 行 式 多 父 代 正 交 交 叉 ， 较 大 程度 上 提 


高 了 算法 的 搜索 效率 。 同 时 ， 为 了 平衡 算法 的 探索 与 开发 能 


力 ， 在 参与 交叉 的 个 体 选 取 与 变异 方程 的 引 


三 


学 习 环节 


别 引 入 动态 自 适应 调节 策略 。 最 后 通过 对 机 


试 与 比较 ， 实 验 结果 表明 ， 改 进 的 算法 
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准 测试 函数 的 测 
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与 以 及 较 高 的 寻 优 精 度 。 
1 ”回溯 搜索 优化 算法 


2 ”元 胞 自动 机 与 正 交 实验 设计 


2.1 元 胞 自动 机 
可 济 搜 索 优化 算法 是 一 种 基于 种 群 迭代 的 智能 进化 算 元 胞 自动 机 055- 《CA) 是 由 Neumann 首先 提出 的 一 种 
法 ， 与 一 般 的 进化 算法 相似 ， 都 上 共有 三 个 基本 的 基因 算 子 ， 经 典 的 自然 计算 模型 ， 具 有 模拟 复杂 系统 时 空 演化 过 程 的 能 
即 交 叉 、 变 异 和 选择 。 有 具体 包括 五 个 步 又 ， 分 别 为 种 群 初始 ” 力 。CA 共 由 元 胞 、 状 态 、 邻 居 和 规则 四 部 分 组 成 ， 是 一 种 
化 、 选 择 T、 种 群 变异 、 种 群 交叉、 选择 I 。 时 间 、 空 间 和 状态 都 具有 离散 性 的 动力 学 系统 。 所 有 的 元 胞 


1) 种 群 初始 化 构成 一 个 元 胞 空间 , 元 胞 空间 中 的 每 个 元 胞 将 并 行 同步 演变 ; 
BSA 的 种 群 初始 化 以 随机 产生 的 方式 进行 ， 且 包括 种 群 包 在 某 一 时 刻 只 能 存在 一 种 状态 ， 且 该 状态 来 源 于 国 
Pop 和 历史 种 群 oldPop 的 同时 初始 化 。 即 ; 眼 离散 状态 集 ， 邻居 是 元 胞 个 体 周围 按 一 定 几何 规则 划分 
Bp 集合 ， 每 个 元 胞 的 下 一 时 刻 状态 由 该 元 胞 以 及 与 其 对 

ee 4 应 的 邻居 共同 决定 ， 规 则 是 根据 元 胞 当前 状态 及 其 邻居 状态 

定 下 一 时 刻 该 元 胞 状态 的 动力 学 函数 。 


n 
加 tH 
三 
于 


屋 冉 侣 对 Hl 
量 


其中 ，ie{1,2.…,SN} ，je{1,2.…,D} SN 和 D 分 别 为 种 群 的 规 


模 和 维 数 ;， low 和 wp 分 别 指 搜索 空间 的 下 界 和 上 界 ;U 为 随 二 维 元 胞 自动 机 的 邻居 模型 划分 通常 可 分 为 四 种 : 冯 - 

机 均匀 分 布 函数 。 诺 依 曼 型 ， 摩 尔 型 ， 扩 展 的 摩尔 型 ， 马 哥 勒 斯 型 。 其 中 最 常 
2) 选择 I 用 的 是 冯 - 诺 依 曼 (Von Neumann ) 型 和 摩尔 (Moore) 型 ， 
BSA 的 选择 工 策略 主要 用 于 确定 历史 种 群 odPop 的 更 以 常用 规则 四 方 网 格 划 分 为 例 ， 如 图 1 所 示 。 


新 。 如 式 (2) 和 (3) 所 示 。 


Pop ifa<b 
et oldPop = . (2) 
和 oldPop otherwise 
” oldPop = randperm(oldPop) (3) 


其 中 : a 和 4b 均 为 服从 区 间 (0,1)〉 上 的 均匀 分 布 随机 数 ; 
randperm 为 重新 排序 函数 。 在 每 次 迭代 过 程 中 ， 历 史 种 群 的 


个 体 都 有 机 会 取 到 当代 以 及 之 前 任意 一 代 的 值 ， 和 凭借 其 独特 (a) Von Neumann 型 (b) Moore 型 
的 记忆 功能 引导 种 群 的 后 续 进化 。 1 常用 元 胞 自动 机 邻居 模型 
3) 种 群 变异 Fig.1 Common cellular automaton neighbor models 
BSA 的 变异 过 程 如 式 (4) 所 示 。 由 于 元 胞 空间 理论 上 在 各 维 上 无 限 延 伸 。 因 此 ， 对 于 二 
Mutant = Pop + F» (oldPop — Pop) (4) ” 维 元 胞 空间 需 定义 周期 型 边界 条 件 , 即 上 下 相 接 , 左右 相 接 ， 
其 中 : Mutant 为 生成 的 变异 种 群 ，F 为 变异 尺度 系数 ， 用 于 ”形成 一 个 拓扑 的 圆 环 面 。 如 图 2 所 示 。 


控制 引导 搜索 向 量 (ol4dPop -Pop) 的 幅度 ，f =3*randn ，randn 为 
服从 标准 正 态 分 布 的 随机 数 。 
4) 种 群 交 叉 
BSA 的 交叉 过 程 包含 两 种 交叉 方式 ,每 个 个 体 通过 对 两 
种 方式 的 随机 调用 ， 生 成 0-1 映射 矩阵 mwaP ， 初 始 wap 为 一 
全 1 和 矩阵 。 有 具体 生成 公式 如 式 (5) (6) 所 示 : 


MAPi nkl : [miralesrand.D) 二 0 (5) 


mapi wana(p) = 0 (6) 
其 中 :w=randperm(D) ，T | 表示 向 上 取 整 ，mixrate 为 设置 的 
交叉 概率 ， rand 为 生成 一 个 [ 叫 区 间 的 随机 数 ; randi(D) 为 


生成 一 个 由 2] 区 间 的 随机 整数 。 四 2 ”二 维 元 胞 空间 拓扑 结 梳 
再 由 确定 的 map 上 映射 矩阵 ， 实 现 种 群 Pop 和 变异 种 群 Fig.2 2-dimensional cellular spatial topology 
Mutant 之 间 的 对 应 交叉 , 最 终 产生 新 的 实验 种 群 7 。 如 式 (7) ”2.2 正 交 实验 设计 
所 示 。 正 交 实验 设计 47 (OED) 是 一 种 利用 正 交 表 来 解决 多 
s a map,.; (7) 姑 素 、 多 水 平实 验 的 有 效 方法 。 在 给 定 的 实验 范围 内 ， 通 过 
9 正 交 性 准则 选取 合适 的 正 交 表 5 (2") ,来 安排 较 少 的 实验 找 


5) 选择 工 
BSA 的 选择 开 策 略 采用 了 贪 焚 选 择机 制 ,在 新 生成 的 实 
验 种 群 7 和 种 群 Pop 之 间 选 择 对 应 位 置 上 的 最 优 (适应 度 值 


最 优 的 水 平 组 合 方式 。L (8" ) 表 示 N 因素，2 水 平 的 正 交 


| 


较 小 的 ) 个 体 用 于 下 一 轮 进化 。 如 式 (8) 所 示 。 表 ，M 表示 需要 安排 的 实验 次 数 。 当 N 和 2 较 大 时 ， 执 行 全 
ks | Ti fimess(T) < fimess(Pop,) (8) 。 面 实验 的 次 数 为 0" , 而 OED 方法 只 需 从 0" 次 全 面 实验 中 抠 
Pop, otherwise 选 出 M 个 均匀 分 散 旦 具有 代表 性 的 点 进行 实验 ( M 通常 远 
重复 以 上 步骤 ， 直 至 满足 循环 终止 条 件 ， 最 后 输出 最 优 小 于 2"” )， 较 大 程度 的 提高 了 工作 效率 。 其 中 ，M =0 ，2 为 

解 。 素数 ，J 为 满足 式 〈9) 的 正 整数 。 
N< = (9) 


QO-1 
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例如 一 个 


的 实验 , 可 以 通过 了 


来 执行 ， 如 式 (10)。 但 若 进 


3 ”改进 的 BSA 
3.1 基于 元 胞 自动 机 的 多 父 代 正 交 交叉 算 子 


> 


T 


Ww = = = 
[> 
一 DD 一 wm PP 


E 交 表 (3’) 


看 实验 则 需 81 (34) 次 。 


(10) 


省 本 回溯 搜索 算法 


又 ， 而 二 项 式 交 叉 算 子 随机 产 台 
应 搜索 空间 的 充分 姑 
实验 设计 的 方法 ， 对 多 个 父 代 进 


的 交叉 方式 是 常用 的 二 项 式 交 
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G, (= (A Py Py ) 
GO)=(P Ba 021. By, ) (14) 


G,(@)= (Bos Por 2. Pon,) 


e) 构 造 正 交 表 Lu (@")=[a, yw ? 根据 正 交 表 Ly (0") 安排 
实验 ， 产 生 W 个 子 代 个 体 〈M 组 实验 ): 


(GO 二 0 Gy (wx)) 
(G (ac (aa Gv(aw)) 


(G(a),cC (an2) Gy (avv)) 


f) 采用 式 (16) 计算 分 析 每 种 因素 的 每 个 水 平 对 实验 结 
果 的 影响 ， 从 而 得 出 最 优 的 水 平 组 合 方式 。 


用 正 交 表 产生 


而 在 确定 的 样本 空间 
模型 ， 提 出 一 种 基于 元 胞 自动 机 
算 子 建立 在 Von Neumann 型 元 胞 
围 的 四 邻居 共同 组 成 多 父 
最 优 的 水 平 组 合 个 体 作 为 该 元 胞 总 
自动 机 的 并 行 演 变 策略 和 


合 元 胞 自动 机 
的 多 父 代 正 交 交 叉 算 子 。 该 


基础 上 ， 每 个 个 体 与 周 


人 参与 正 交 交叉 ， 最 后 


子 代 中 


。 由 于 元 胞 
生 能 ， 使 得 该 


能 性 ， 在 交叉 算 子 中 引入 了 随机 扰动 策 
交叉 概率 时 执行 多 父 代 了 
选任 意 维 以 初始 化 


样 性 。 


新 型 交叉 算 子 在 寻 优 速度 外 
同时 ， 为 了 


后 期 陷入 局 部 最 优 而 


于 元 


9 较 大 的 提升 。 
EE 熟 收 合 的 可 


E 交 交叉 ， 反之 ， 则 对 
的 方式 进行 扰动 。 进 


算 子 操作 


a) 判断 元 胞 个 体 是 否 满足 交叉 概率 mixrare 。 
rand < mixrate 成 立 则 跳 转 至 b)， 执 行 多 父 代 了 


对 该 个 体 进行 随机 扰动 : 


有 1 +( 


b) 将 参与 交叉 的 2 个 父 代 个 体 分 别 作为 了 
一 个 水 平 ，。 记 第 7 外 水 纪 


x Ee {1 为 a xo} 


c) 当 优 化 问题 的 维 


up; —low, )*rand if j erandperm(D,dim) 


元 胞 个 体 满足 
该 个 体 随 机 挑 
虽 了 种 群 的 多 


如 果 


E 交 交叉 。 和 否则 


(11) 
(12) 


E 交 实验 设计 的 
则 B=%, 


匹配 ， 但 是 通常 在 复杂 的 高 维 约束 优 


l1<k<k<...<ke<D, 


则 需 对 维度 进行 切片 划分 ， 从 而 降低 OED 的 了 
将 参与 交叉 的 任 一 个 体 半 分 成 ECE=N) 个 片段 ， 每 个 片段 即 


对 应 一 个 因素 ， 如 式 (13〉 所 示 。 


中 ,kKe 为 (1,D) 之 间 随机 入 


维度 进行 
素 的 选取 
[ 作 量 。 假 定 


(13) 


满足 


d) 令 后 =0, 丰 =N ,因此 第 j (07=12…N) 个 因素 的 2 个 水 


平 可 以 表示 为 


ps 
S$, = Ey (16) 
不 能 实现 对 相 
E 交 交叉 算 子 48-20 利 用 正 交 其 中 : 丈 为 第 mm=L12…M) 组 实验 的 结果 ，S% 为 因素 
的 划分 ， 并 采 za(2=12.…N) 的 第 gq(q=1,2,.…,0) 个 水 平 对 结果 的 影响 程度 。 
E 成 子 代 个 体 ， 从 zm 表示 第 mm 组 实验 中 是 否 包 含 因素 7 的 第 4 个 水 平 , 若 包 含 


则 zm -1， 否 则 zw 0。 
3.2 交叉 选择 机 制 

BSA 算法 的 交叉 环节 中 ， 通 过 map 映射 矩阵 完成 对 变异 
种 群 Muar 和 种 群 Pop 的 随机 交叉 操作 。 昌 然 已 将 原始 交叉 
葵 换 为 基于 元 胞 自动 机 的 多 父 代 正 交 交叉 算 子 ， 较 好 地 提高 
了 算法 的 搜索 效率 和 精度 ， 但 在 搜索 方向 上 仍然 存在 一 定 的 
言 目 性 。 为 此 ， 受 遗传 算法 交叉 前 进行 一 定 概率 的 择优 选取 
启发 ii， 本 文 引 入 交叉 选择 机 制 ， 并 提出 一 种 动态 优秀 个 体 
比例 权重 p% 。 用 于 对 参与 交叉 的 个 体 进行 选择 ， 在 算法 执 
行 的 不 同 阶段 以 合适 的 搜索 方向 引导 交叉 。 
当 种 群 变异 结束 后 ， 变 异种 群 Muar 和 种 群 Pop 将 作为 
候选 交叉 成 员 。 从 候选 成 员 中 选取 适应 度 值 最 优 的 前 Swvp% 
个 个 体 ， 再 从 剩 下 的 成 员 中 随机 选取 sw(I- pg%) 个 个 体 共同 
组 成 交叉 群体 。 优 秀 个 体 比例 权重 p% 满足 如 下 公式 : 

p= pant (pra Pon)" 2S] 47 

中 : ps 和 po 分别 为 权重 的 最 小 值 和 最 大 值 。 epk 和 epoch 
分 别 为 算法 的 当前 适 代 次 数 和 选 代 总 次 数 。 由 式 (17) 可 以 
看 出 ， 随 着 进化 代数 的 不 断 增 加 ， 优 秀 个 体 比例 权重 p% 起 
来 越 大 。 使 得 算法 在 执行 前 期 交叉 算 子 偏向 于 全 局 搜索 ， 而 
到 了 中 后 期 时 则 偏向 于 在 全 局 最 优 值 附近 进行 精确 的 局 部 搜 
索 。 在 搜索 方向 上 起 到 了 较 好 的 引导 作用 。 
3.3 动态 变异 方程 
通过 对 BSA 算法 的 变异 算 子 分 析 可 知 , 该 算法 的 变异 操 
作 能 记忆 之 前 迁 代 过 的 种 群 ， 并 以 迭代 过 的 历史 种 群 作为 经 
验 来 指导 当前 种 群 的 搜索 方向 ， 使 之 具有 较 强 的 全 局 搜索 能 
力 。 但 因 其 仅仅 通过 设置 历史 种 群 对 搜索 方向 进行 引导 ， 个 
体 进化 缺乏 一 定 的 学 习 能 力 ， 导 致 算法 的 收敛 速度 减 慢 ， 局 
部 开采 能 力也 较 弱 。 针对 这 一 缺点 ,在 文献 [5.13] 提 出 的 广泛 
学 习 策略 以 及 ABC 算法 雇佣 峰 搜索 方程 的 基础 上 ， 本 文 
在 BSA 的 变异 过 程 中 同时 加 入 最 优 个 体 和 较 优 群体 的 引导 ， 
提出 一 种 带 有 加 权 因 子 的 自 适应 动态 变异 方程 。 公 式 如 式 
(18) (19) 所 示 。 


Mutant, =u ( Pop;+F “(oldPop, — Pop,, )) 十 


(1—u)(Popys, + FoldPop,, ww 一 Pop)) (18) 
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hina 
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录用 


定稿 


w=exp(~ape-1) 


其 中 : 4 为 加 权 因 子 , i, ,为 [hbSN] 上 互 不 相等 的 随机 整数 ; 
局 最 优 个 体 ，ol4dPop,, we 为 当代 
历史 种 群 中 前 +(mod(SN,7)=0) 个 优秀 个 体 组 成 的 较 优 群体 。 


Popu 为 当前 种 群 搜索 到 的 全 


选 出 的 种 群 ol4Pop, js 和 种 群 0 


杨 军 ， 等 : 


(19) 


ldPop 规模 相同 ， 具 体操 作 如 


下 : 将 oldPop 中 选 出 的 前 1 个 优秀 个 体 复制 成 SNVt 组 ， 并 对 


每 一 组 随机 排序 ， 最 后 重新 组 合成 种 群 ulPopv ur 。 
式 (18) 的 前 半 部 分 通过 选取 种 群 odPop 和 Pop 中 与 变 
差分 向 量 来 增强 算 子 的 勘探 能 力 ， 


异 个 体 i 不 同 的 个 体 组 成 差分 
后 半 部 分 通过 全 局 最 优 个 
oldPop。,scs 来 增强 算 子 的 精英 学 
式 (19) 中 在 算法 执行 的 不 同 


力 进行 合适 的 加 权 。 使 得 算法 执行 前 


后 期 则 以 局 部 开发 为 主 。 
3.4 改进 算法 流程 


阶段 对 全 局 搜索 和 


一 种 基于 元 胞 自动 机 的 动态 回溯 搜索 优化 算法 

与 交叉 个 体 的 选择 
媳 将 选 旧 
axa(axa=SN) 的 二 维 网 格 元 胞 空间 中 。 


的 
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参与 交叉 群体 随机 分 布 在 


规模 为 


多 执行 交叉 操作 : 以 交叉 概率 maoare 对 元 胞 空间 中 的 每 


体 Popss 和 较 优 历史 群体 
习 能 力 和 开采 能 力 。 并 且 利 用 
局 部 搜索 能 


局 搜索 为 主 ， 而 


期 以 全 


epk = epk 


进化 的 优秀 种 群 Pop 。 


ji 找 出 第 epk 代 的 最 优 个 体 和 最 优 适 应 度 值 。 


j) 淹 


区 输出 当前 的 最 优 个 体 和 最 优 适应 度 值 。 


个 元 胞 个 体 进行 判断 。 如 果 满 足 则 以 Von Neumann 型 元 胞 邻 
居 结 构 为 基础 ， 构 造 正 交 表 (8" ) 执 行 多 父 代 正 交 交 又 ,并 
获得 最 优 组 合 形式 的 子 代 。 否 则 以 式 (11) (12) 对 个 体 进行 
随机 扰动 。 最 终 得 出 实验 种 群 7 。 

I) 执 行 选 择 开 策略 : 评价 实验 种 群 T ， 选 出 用 于 下 一 轮 


断 是 否 达 到 循环 终止 条 件 。 若 是 则 转 步骤 七， 否则 
+1， 跳 转 至 步骤 b)。 


a 8 4 ”数值 实验 及 结 
改进 后 的 CA-DBSA 算法 具体 实现 流程 如 下 : 数值 实验 及 结果 分 析 
a) 初 始 化 种 群 Ppp 、 历 史 种 群 o14Pop 以 及 种 群 规模 SN 、 为 了 验证 本 文 算法 的 可 行 性 与 效果 ,选取 12 个 标准 测试 
交叉 概率 mixrate 和 最 大 进化 代数 epoch 等 参数 , 设置 当前 进化 函数 [565 对 其 进行 数值 实验 。 同 时 与 基本 BSA 算法 QI、 参考 
代数 epk=1 。 文献 中 的 CLPSO51、MABCBI 和 LBSA04 算 法 ， 以 及 利用 正 
b) 执 行 选择 I 策略 : 更 新 历史 种 群 oldPop 。 交 交 又 策略 改进 的 OXDE 和 多 父 代 交叉 算 子 的 
c) 评 价 种 群 Pop 和 历史 种 群 o14Pop 的 适应 度 值 , 选 出 最 优 MPCO-ABCDB9 算 法 进行 性 能 比较 。 
个 体 Pom 和 较 优 历史 群体 "ldPopw ww 。 测试 函数 如 表 1 所 示 。 表 中 给 出 了 函数 具体 名 称 和 表达 
dg) 执行 变 异 操作 : 根据 式 (18) (19) 计算 出 加 权 因 子 4 式 、 搜 索 空间 范围 以 及 最 优 值 。 其 中 ~ 三 是 单 峰 类 型 函数 ， 
和 变异 种 群 Mutant 。 侧重 于 衡量 算法 的 收敛 速度 和 精度 ; 有 ~ fs 是 多 峰 类 型 函数 ， 
e) 执 行 交 叉 选 择机 制 : 根据 式 (17)〉 计算 出 当代 优秀 个 侧重 于 衡量 算法 的 全 局 搜索 能 力 。 
体 比 例 权 重 P% ,并 从 种 群 Pop 和 变异 种 群 Mutant 中 完成 对 参 
表 1 12 个 标准 测试 函数 
‘Table 1 12 standard test functions 
函数 名 函数 表达 式 搜索 空间 最 优 值 
Sphere fi(x) 3 [-100,100] 0 
Tablet f(x)=10°0 + [-100,100] 0 
Rosenbrock f(x)= S(100(%, a ) +(1—% )) [-30,30] 0 
Schwefel 1.2 f(x) yb [-100,100] 0 
Schwefel 2.22 fs(7) =>hl+ TI [-10,10] 0 
Ackley f(x)=20+e— exp [02 E> | 一 exp (Deos(rs )] [-32,32] 0 
Alpine 万 (加 => 所 sin(%)+0.1| [-10,10] 0 
Rastrigin OPSAE —10cos(2xx,)+10) [-5.12,5.12] 0 
Griewank f(x) -六 -TIee( 却 j [-600,600] 0 
Schwefel 2.26 fi (2) =418.98288727216249D — > (% sin( Vs|)) [-500,500] 0 
fii(x)= (os (oy) + -1) (1 +10sin (xy, +1)) + 
Levy and Montalvo 1 n 1 [-10,10] 0 
(39,—1)),», =1+0.25(x, +1) 
J (x)= 0.1(sin? (3xx,) + (x -1) (1 + sin’” (3xx, )) + 
Levy and Montalvo 2 四 [-5,5] 0 


(xx 一切 (1 +sin’ (2zx, ))) 


本 次 实验 中 , 7 种 对 


mixrate 设 为 1; 


别 为 10、90; 正 交 表 采 用 六 


杨 军 ， 


等 : 一 种 基于 元 胞 自动 机 的 动态 回潮 


比 算法 的 种 群 规模 SN 和 最 大 i 
数 epoch 均 设 置 为 25 和 5000; 基本 BSA 算法 的 交叉 概率 
CA-DBSA 中 优秀 个 体 比 例 权 重 


因素 五 水 平 的 (5°) 。 其 


Pnmin 、 


ore 7,9,14,20] 进 行 对 应 设置 。 


Q=0.05 的 Wilcoxon 秩 和 检验 方法 对 其 结果 六 
-、 罕 分 别 表 示 
CA-DBSA 算法 。 


实验 分 为 三 


行 , 第 1 


二 


其 他 算法 的 诉 


E 能 要 优 于 、 


、2 组 实验 中 , 测试 函数 自 
D 均 取 30; 第 3 D 均 取 60。 针 对 
仿真 实验 , 每 种 算法 独立 
均值 (Mean) 和 标准 差 (Std.Dev)。 为 了 进 
种 算法 之 间 的 性 能 差异 , 在 第 


化 代 
况 来 看 ， 当 


搜索 优化 算法 


由 算法 


Py 


身 维 数 
同一 测试 函数 的 
运行 30 次 , 记录 其 寻 得 最 优 值 的 平 


步 分 析 
2、3 组 实验 中 引入 显著 1 


行 直观 呈 


劣 于 、 相 


较 7 


生 水 平 


现 。+、 


当 于 


寻 优 能 力 在 12 个 测试 函 
时 ， 算 法 开采 能 力 较 弱 ， 


则 易 陷 入 局 部 最 全 
本 文 的 CA-DBSA 
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优 测试 结果 已 对 其 进行 加 粗 体现 。 从 最 优 实 验 结果 的 分 布 情 
当 交 叉 概 率 mixrate 取 值 为 0.7 时 , CA-DBSA 算法 的 


数 上 均 表现 最 好 。 当 交叉 概率 偏 小 


算 法 交 又 概 率 mixrate 


表 3 给 


优 结果 是 


CA-DBSA 算法 在 大 部 


他 几 种 算法 ， 表 现 出 了 较 强 的 全 局 搜索 和 局 部 搜索 能 力 。 相 
的 CLPSO、MABC 和 


对 于 改 i 


局 部 搜索 不 充分 。 反 之 ， 在 偏 
ti。 因此, 在 接 下 来 的 第 二 和 第 三 组 实验 中 ， 


大 时 


UD 


Li y ~ 


均 取 值 为 0.7。 
出 了 第 2 组 实验 D=30 时 各 算法 的 测试 结果 。 从 
表 中 每 个 函数 的 最 优 测试 结果 可 以 看 出 
CA-DBSA 算法 实验 得 到 。 表 明了 本 文 的 
分 测试 函 数 上 的 实验 结果 都 要 优 于 其 


有 9 个 函数 的 最 


基本 BSA 来 说 ， 


CA-DBSA 


算法 在 寻 优 精度 和 稳定 性 上 均 具 有 明显 的 优势 ; MPCO-ABC 


和 LBSA 均 分 别 只 在 一 个 函数 上 性 能 优 于 CA-DBSA; 而 对 


为 了 证 明 本 文 所 提 算 法 CA-DBSA 交叉 概率 mixrate 的 最 。 于 OXDE， 只 有 在 函数 请 、 大 和 访 上 ，OXDE 算法 的 寻 优 
理想 取 值 ， 在 第 一 组 实验 中 ， 对 mixrate 分 别 取 0.5~0.9 的 5 能 力 比 本 文 算法 更 优 ， 且 上 相当 于 本 文 算法 。 这 正 也 体现 
个 不 同 的 值 进行 测试 。 结 果 如 表 2 所 示 ， 其 中 每 个 函数 的 最 。 出 了 正 交 实 验 设 计 方法 的 有 效 性 。 

表 2 时 ca-dbsa 算法 mixrate 不 同 取 值 的 实验 结果 
Table 2 CA-DBSA algorithm with D=30 experimental results with different values of mixrate 
Mean+t Std.Dev 
函数 mixrate =0.5 mixrate =0.6 mixrate =0.7 mixrate =0.8 mixrate =0.9 
2.82E-114 十 3.41E-114 7.03E-118 土 3.75E-117 3.56E-121 士 1.38E-120 2.29E-120 土 3.42E-119 1.63E-116 土 5.97E-117 
太 1.17E-112 圭 5.26E-110 1.85E-117 十 2.44E-116 1.60E-122 土 4.73E-122 3.04E-119 士 3.26E-119 6.05E-115 士 3.09E-116 
A 7.44E+01 十 5.32E+01 5.81E+01 士 2.47E+01 1.81E+00 土 1.56E-01 2.23E+00 土 1.07E+00 3.96E+01 土 1.20E+01 
万 5.41E-158 土 3.11E-160 2.63E-165 圭 2.79E-164 1.47E-174 土 3.10E-174 1.83E-171 士 4.36E-171 1.47E-164 士 4.72E-166 
fs 3.02E-63 土 1.48E-63 4.92E-68 土 1.89E-66 1.37E-75 土 8.77E-75 5.26E-72 土 7.25E-71 2.32E-66 二 6.74E-66 
J 9.63E-09 士 9.91E-08 5.03E-11 土 4.99E-10 2.03E-12 士 4.S0E-12 4.15E-10 十 3.89E-10 3.37E-09 土 4.68E-10 
万 4.81E-11 土 6.23E-12 3.47E-12 土 1.16E-12 7.94E-15 土 1.03E-15 5.11E-14 土 1.57E-13 2.87E-12 土 4.1SE-11 
kk 6.05E-01 土 4.67E+00 4.76E-01 土 2.67E-01 2.05SE-01 士 3.42 卫 -02 8.49E-01 士 6.64E-01 5.29E+00 土 1.10E+00 
h 1.93E-07 十 2.36E-06 0.00E+00 土 0.00E+00 0.00E+00 土 0.00E+00 2.05E-10 土 1.35E-11 2.29E-08 土 4.51E-08 
fio 5.83E-09 士 2.26E-10 1.52E-11 士 1.47E-11 1.82E-12+1.13E-13 9.01E-10 士 7.23E-11 7.05E-10 土 4.43E-10 
fi 2.76E-16 圭 3.71E-16 4.73E-19 土 1.05E-20 1.49E-21++1.75E-22 1.19E-18 士 1.32E-18 3.54E-17 土 2.16E-16 
fa 6.87E-24 圭 3.61E-26 3.13E-27 土 3.65E-28 1.13E-31 土 2.68E-31 3.24E-29 土 4.41E-29 1.17E-25 十 2.09E-26 
表 3 时 7 种 算法 的 实验 结 
Table 3 Experimental results of 7 algorithms at D=30 
函 数 CLPSO MABC MPCO-ABC OXDE BSA LBSA CA-DBSA 
Mean 2.43E+00 1.42E-17 3.69E-16 2.02E-26 9.36E-20 5.88E-71 3.56E-121 
. Std.Dev 1.92E-01 6.89E-17 5.68E-15 1.23E-25 4.64E-20 3.32E-68 1.38E-120 
Mean 5.90E+00 1.13E-18 1.58E-16 2.79E-27 2.51E-18 4.45E-73 1.60E-122 
Std.Dev 8.46E+00 4.84E-18 2.91E-14 9.11E-26 2.24E-18 2.24E-72 4.73E-122 
Mean 1.24E+02 2.88E+01 3.60E-03 1.11E-02 2.50E+01 1.30E-01 1.81E+00 
’ Std.Dev 5.23E+02 1.27E+00 8.03E-03 1.34E-02 3.21E+01 2.52E-01 1.56E-01 
Mean 1.09E+04 6.89E-14 5.09E+03 1.58E-48 4.00E-36 4.13E-139 1.47E-174 
有 Std.Dev 6.04E+03 2.37E-15 3.81E+03 4.65E-48 4.73E-35 5.68E-137 3.10E-174 
Mean 1.79E-01 2.40E-10 1.20E-15 7.75E-15 3.14E-14 2.12E-39 1.37E-75 
” Std.Dev 3.47E-01 6.22E-11 7.34E-15 1.26E-16 1.52E-12 1.74E-39 8.77E-75 
Mean 1.95E-01 1.35E-03 3.06E-09 7.19E-14 1.65E+00 1.34E+00 2.03E-12 
* Std.Dev 1.61E-01 1.33E-03 1.26E-10 5.77E-15 1.34E+00 2.86E-01 4.50E-12 
Mean 5.38E+00 3.80E-05 2.69E-08 3.24E-12 7.44E-14 1.92E-12 7.94E-15 
7 Std.Dev 2.18E-01 2.54E-07 2.12E-08 4.65E-13 1.42E-14 4.78E-12 1.03E-15 
Mean 3.21E+02 9.02E+01 6.29E-01 9.95E+00 1.39E+01 1.03E+01 2.05E-01 
A Std.Dev 9.16E+01 1.50E+00 1.07E+01 6.03E-01 2.19E+01 1.48E+00 3.42E-02 
Mean 6.75E-01 4.92E-07 2.01E-16 0.00E+00 1.11E-16 9.90E-03 0.00E+00 
» Std.Dev 1.28E-02 9.16E-08 4.43E-16 0.00E+00 4.47E-16 6.45E-01 0.00E+00 
fio Mean 8.02E+01 4.32E-04 7.59E-06 1.18E+02 1.66E+03 1.95E+03 1.82E-12 
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Std.Dev 1.52E+01 4.09E-03 


Mean 1.04E+01 4.95E-01 
Std.Dev 1.12E+01 1.60E+00 
Mean 1.36E+00 3.60E-05 


Std.Dev 1.38E+00 2.09E-04 
+ 0 0 
- 12 12 
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除了 实验 结果 的 寻 优 精度 和 稳定 性 ， 收 敛 速度 也 是 衡量 
算法 能 力 的 一 项 重要 指标 。 图 
个 测试 函数 上 优化 过 程 的 收敛 


过 程 


日 


3 给 出 了 D=30 


> 


线 的 变化 ， 图 中 的 纵 


TT 


线 。 为 了 更 清晰 地 描述 收敛 
标 取 适 应 度 值 的 对 数 。 从 图 
可 以 看 出 ， 在 没有 早熟 收敛 的 前 提 下 ， 与 其 他 6 种 算法 相 
比 ，CA-DBSA 算法 在 12 个 函数 上 整体 来 说 均 表 现 出 了 较 快 


闻 


时 各 算法 在 


jh 
[We 


的 收敛 速度 。 除 了 在 函数 大 上 ， 本 文 算 法 迭代 到 5000 代 时 
仍 未 完全 收敛 ， 在 其 他 函数 上 均 在 前 半 周 期 已 经 收敛 。 这 得 


TI 


提高 。 


益 于 基于 元 胞 
CA-DBSA 在 搜索 效率 上 相 比 于 基本 BSA 得 到 了 较 大 程度 的 
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对 于 智能 进化 算法 而 言 , 随 着 优化 问题 维度 的 不 断 增 加 ， 从 表 中 加 粗 部 分 数据 结果 和 秩 和 检验 的 统计 信息 来 看 ，12 
问题 的 复杂 度 也 将 会 随 之 急剧 上 升 。 而 对 于 多 峰 类 型 函数 来 ”个 函数 的 最 优 测试 结果 中 有 10 个 是 由 CA-DBSA 算法 实验 所 
说 ， 其 局 部 极 值 个 数 将 呈现 指数 级 别 的 增长 。 因 此 算法 的 整 得 ， 比 表 3 实验 结果 表现 更 好 。 在 OXDE 中 ,只 有 在 函数 
本 性 能 将 会 受到 较 大 影响 。 在 此 基础 上 ， 为 了 进一步 研究 本 和 上 ，OXDE 的 寻 优 结果 比 CA-DBSA 更 好 。 而 表 3 中 存 
文 所 提 算 法 受 维度 增长 的 影响 程度 ,对 7 种 算法 在 12 个 测试 “在 4 个 结果 比 CA-DBSA 更 优 或 相似 。 同 样 ,30 维 时 的 LBSA 
函数 上 进行 了 第 三 组 D=60 时 的 数值 实验 。 实 验 结果 如 表 4 “有 一 个 函数 寻 优 性 能 优 于 CA-DBSA， 而 在 60 维 时 性 能 有 所 
所 示 。 从 表 中 可 以 看 出 ， 当 测试 函数 的 维 数 从 30 增加 到 60 ”下降 。 因 此, 与 表 3 对 比 可 发 现 ， 本文 算法 在 60 维度 时 的 整 


时 ， 相 比 于 表 3 中 的 实验 结果 ， 和 名 算法 的 收敛 性 能 均 有 所 下 体 实 验 结果 比 30 维度 时 有 所 提高 ， 展 现 出 了 较 强 的 鲁 棒 性 。 
降 , 但 是 CA-DBSA 受到 的 影响 相对 其 他 6 种 算法 来 说 较 小 。 
表 4 时 7 种 算法 的 实验 结果 
Table 4 Experimental results of 7 algorithms at D=60 
函数 CLPSO MABC MPCO-ABC OXDE BSA LBSA CA-DBSA 
Mean 1.33E+02 7.27E-07 5.51E-09 7.07E-11 6.63E-06 2.07E-24 1.33E-52 
有 Std.Dev 4.17E+01 9.89E-08 4.62E-08 7.90E-11 3.90E-05 1.52E-26 7.29E-52 
Mean 1.46E+02 2.37E-06 9.04E-08 5.27E-10 7.01E-06 5.96E-27 $.77E-S2 
有 Std.Dev 2.25E+00 5.76E-05 6.61E-08 1.19E-11 9.96E-06 3.62E-28 1.20E-53 
Mean 3.13E+04 5.73E+04 7.48E+02 1.56E+02 2.41E+02 4.02E+02 8.73E+01 
Std.Dev 1.94E+03 1.41E+04 4.23E+01 2.54E+02 2.23E+02 3.68E+01 3.48E+00 
Mean 1.44E+06 4.29E+02 5.16E+05 2.91E-19 2.40E-09 7.62E-31 5.81E-69 
4 Std.Dev 6.59E+03 1.28E+02 6.49E+04 5.64E-20 3.47E-09 1.26E-28 1.60E-70 
Mean 5.02E+03 2.72E-04 7.22E-07 4.94E-05 8.79E-04 3.79E-08 4.28E-35 
万 Std.Dev 7.74E+01 2.21E-06 6.36E-08 8.80E-05 6.50E-03 1.17E-06 1.1SE-34 
Mean 2.04E+01 1.52E+01 2.67E-03 2.52E.-06 5.58E+00 6.40E+02 2.18E-02 
4 Std.Dev 2.01E+01 8.05E+00 1.83E-04 1.S0E-0S 3.25E-02 3.31E+02 2.66E-04 
Mean 1.08E+02 6.43E-01 8.04E-05 9.78E-07 1.19E-04 8.14E-09 1.83E-14 
1 Std.Dev 1.05E+02 6.78E-01 1.62E-03 1.52E-08 3.76E-05 6.22E-08 1.52E-16 
Mean 7.87E+02 1.75E+02 4.37E+01 5.17E+01 4.48E+01 8.94E+01 2.07E+01 
有 Std.Dev 1.03E+01 7.58E+02 2.03E+00 6.17E+01 4.74E+01 6.25E+00 2.56E+00 
Mean 8.72E+02 2.47E-03 2.41E-06 3.11E-11 9.87E-03 1.30E+01 5.13E-10 
’ Std.Dev 3.07E+02 5.95E-04 3.72E-06 3.36E-10 1.62E-01 5.34E+01 7.12E-10 
Mean 1.63E+04 1.05E+04 7.81E+02 3.77E+03 3.57E+03 7.71E+03 1.37E+02 
加 Std.Dev 1.87E+04 8.87E+03 6.33E+02 3.57E+03 3.95E+03 5.22E+03 2.37E+00 
Mean 1.66E+02 1.85E+00 5.18E+02 7.10E-02 2.19E-01 1.06E+01 1.2SE-0S 
0 Std.Dev 1.71E+01 7.51E-01 1.58E+01 9.98E-02 2.23E-01 2.05E+00 4.69E-04 
Mean 4.14E+01 5.89E-01 1.73E-09 6.33E-14 3.28E-08 1.49E+01 8.71E-18 
如 Std.Dev 3.60E+01 3.23E-02 3.57E-07 1.26E-13 2.10E-07 6.37E+02 7.11E-18 
十 0 0 1 2 0 0 
- 12 12 11 10 12 12 
0 0 0 0 0 0 
5 结束 语 [M].Cambridge: MIT Press, ee | 
[2] Kennedy J, Eberhart R. Particle Swarm optimization [C]//Proc of IEEE 
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